Collective excitations in ferrimagnetic Heisenberg ladders 
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We study ground-state properties and the low-lying excitations of Heisenberg spin ladders com- 
posed of two ferrimagnetic chains with alternating site spins (5*1 > 5*2) by using the bosonic Dyson- 
Maleev formalism and Lanczos numerical techniques. The emphasis is on properties of the ferrimag- 
netic phase which is stable for antiferromagnetic interchain couplings J± > 0. There are two basic 
implications of the underlying lattice structure: (i) the spin-wave excitations form folded acoustic 
and optical branches in the extended Brillouin zone and (ii) the ground state parameters (such as the 
on-site magnetizations and spin-stiffness constant) show a crossover behavior in the weak-coupling 
region < J_l < 1- The above peculiarities of the ladder ferrimagnetic state are studied up to second 
order in the quasiparticle interaction and by a numerical diagonalization of ladders containing up to 
A*' = 12 rungs. The presented results for the ground-state parameters and the excitation spectrum 
can be used in studies on the low-temperature thermodynamics of ferrimagnetic ladders. 



PACS: 75.10.-b, 75.10.Jm, 75.40.Gb, 75.50.Gg 



I. INTRODUCTION 



There is an impressive difference between the proper- 
ties of a single spin-i antiferromagnetic Heisenberg chain 
and the related ladder system composed of two coupled 
chains. The spin-^ chain is critical, with a power-law 
decay of the spin correlations, whereas the spin-i ladder 
has a finite energy gap in the excitation spectrum and ex- 
ponentially decaying correlations .EJ Similar behavior has 
been found in multi-leg spin-i ladders with odd and even 
number of chains, respectively^ Since in spin ladders the 
excitation spectrum is controlled by two energy scales 
(the intrachain exchange constant J and the transverse 
interchain exchange coupling J±), variation of J± and 
the applied magnetic field may produce a rich variety of 
specific quantum effects and phase transitions. In par- 
ticular, intermediate plateaus in the magnetization pro- 
cesses of a number of ladder systems have been predicted 
and experimentally observed in the last few years.cfQ 

Synthesised quasi-one-dimensional mixed-spin 
compoundaa constitute an appropriate base for future 
developments in the physics of spin ladders. Most of 
them are molecular magnets containing two different 
transition-metal magnetic ions alternatively distributed 
on the lattice. Published experimental work implies that 
the magnetic properties of these mixed-spin materials are 
basically described by quantum Heisenberg spin models 
with antiferromagnetically coupled nearest-neighbor lo- 
calized spins. Some typical examples of bipartite ladder 
structures composed of two different spins are shown 
in Fig. ^ The first two ladder structures reproduce, 
e.g., arrangements of the magnetic atoms Mn (5*1 = |) 
and Cu (^2 — 5) along the o-axis in the compounds 
MnCu(pbaOH)(H20)3 (pbaOH = 2-hydroxy-l,3- propy- 
lenebisoxamato) and MnCu(pba)(H20)3-2H20 (pba = 
1,3-propylenebisoxamato), respectively. Along the c- 



axis the magnetic ions in both mixed-spin compounds 
are arranged as shown in Fig. 0(c) Jj Since the molec- 
ular chemistry admits a relatively easy control of the 
molecular-unit positions, it may be expected that the 
discussed mixed-spin ladders will be synthesized in the 
near future. 
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FIG. 1. Typical bipartite mixed-spin ladders composed of 
two types of site spins Si(n) and 82(^1): Si(n) = 51(5*1 -|- 1), 
Si(n) = 52(52 + 1) (5i > 52). 

Ground-state properties of the Heisenberg spin mod- 
els defined on the ladders shown in Fig. |l| have recently 
been studied by using semiclassical nonlinear-, a model 
techniquesB'Q and strong-coupling expansions B'H It was 
found that an interplay between bond and spin alter- 
nations can produce various gapful phases separated by 
critical lines on the phase diagram. In this article we 
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address the mixed-spin ladder model which is shown in 
Fig. [|(a) and defined by the Heisenberg Hamiltonian 



N 



7^ = Si(n) • S2(n + 1) + S2(n) -81(71 + 1) 

n=l 
N 

+ Ji^Si(n)-S2(n). (1) 



The integers n label N rungs, each of them contain- 
ing two different site spins with quantum spin num- 
bers .S*! > 82- Below we frequently use the notations 
CT = S'i/S'2 > 1 and ^2 = S, and set J — 1. 

For antiferromagnetic interchain couplings (Jx > 0), 
the model (^ exhibits a ferrimagnetic ground state chafc. 
acterized by the net ferromagnetic moment (6*1 — S'2)A^.t2l 
Such a magnetic phase may be referred to as a quan- 
tized unsaturated ferromagnetic phase: it is character- 
ized both by the quantized ferromagnetic order param- 
eter M — X]n=i [Si (»^) + 82(71)] [quantized in integral 
or half-integral multiples of the number of rungs, M — 
{Si — S2)N] and by the macroscopic sublattice magne- 

tizations Ma = En=i{S!in)) and Mb = 
In the limit J± ^ 00 the rung spins form local spin 
states with a total spin Stot = Si — 82'- in particular, the 
(1, i) ladder is equivalent to the spin-i ferromagnetic 
Heisenberg model with an effective exchange interaction 
Jeff = — |. On the other hand, for ferromagnetic inter- 
chain couplings {Jj_ < 0), it may be generally expected a 
magnetically disordered spin-singlet ground state, since 
the isotropic Heisenberg model (|l|) is defined on a bipar- 
tite one-dimensional (ID) lattice. 

In spite of the fact that the ferrimagnetic long-range or- 
der exists already in the limit of nonintcracting chains, in 
the antiferromagnetic region {J± > 0) the checkerboard 
ladder exhibits various specific properties deserving a 
special study. As shown below, the most interesting re- 
gion corresponds to relatively small interchain couplings 
{0 < J± < 1) where a number of ground-state parameters 
exhibit a crossover behavior: the latter is marked by the 
extreme values of these parameters. For instance, in that 
region the sublattice magnetization Ma and the spin- 
stiffness constant ps reach their maximal values. An- 
other specific property of the ferrimagnetic ladder con- 
cerns the collective spin-wave excitations, which form in 
the extended Brillouin zone —tt/qq < k < ir/ao folded 
acoustic and optical branches characterized by a mini- 
mum at the zone boundary k = n/ao and a maximum 
between k — 7r/(2ao) (for J± = 0) and k = Tr/ag (for 
J_L =00), ao being the lattice spacing along the ladder. 
A purpose of the present work is to analyze the above 
mentioned peculiarities of the ladder system. 

The existence of a magnetic ground state opens the 
possibility of using systematic spin-wave approaches 
for the ID Hamiltopian (|l|). Indeed, recently pub- 
lished calculationsEj1i3 argued that the spin-wave the- 
ory (SWT) is capable to produce precise quantitative 



results both for the ground-state parameters and exci- 
tation spectrum of ferrimagnetic Heisenberg chains, in- 
cluding the extreme quantum system (1,|). Turning 
to the ferrimagnetic ladders, a SWT description in the 
weak-coupling limit Jj^ <C 1 looks to some extend puz- 
zling since in that region one may expect strong fluctu- 
ations of the macroscopic chain magnetizations around 
the common quantization axis. Similar problems arise 
when the semiclassical nonlinear a modpL technique is 
applied to antiferromagnetic spin ladders .liJ A successful 
treatment of the weak-coupling limit Jj^ <C 1 in ferri- 
magnetic ladders requires an appropriate choice of the 
free-quasiparticle bosonic Hamiltonian. Using as a guide 
the Lanczos numerical technique, it is argued below that 
the Dyson-Maleev formalism is capable to produce pre- 
cise quantitative description of the ferrimagnetic state in 
a relatively large range of interchain couplings, including 
the weak-coupling limit Jj_ <C 1. 

The paper is organized as follows. In Section 2 we in- 
troduce the bosonic representation of Ti and discuss the 
choice of an appropriate free-quasiparticle bosonic Hamil- 
tonian. As a criterion we use the perturbation series for 
Ma compared to Lanczos numerical results. In Section 
3 we study the excitation spectrum of the model up to 
second order in the bosonic interaction V. Section 4 con- 
tains a summary of the results and discussions concerning 
also the disordered phase. In the Appendix we calculate 
the second-order 0{V^) corrections for Ma- 



II. BOSONIC REPRESENTATION OF H: 
PERTURBATION SERIES FOR Ma 

A. Bosonic Hamiltonian 

We adopt the Dyson-Maleev formalism.El Performing 
subsequently the Dyson-Maleev, Fourier, and Bogoliubov 
transformations, we obtain the following Hamiltonian 
terms of the quasiparticle bosonic operators ak and Pk 



nB = Eo+Ho + V, V = V2 + Vdm^ 



(2) 



Eq is the ground-state energy of the ferrimagnetic state 
calculated up to first order in 1/S (see Fig. g): 
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2^ J±/2 + 7fc 

CT+1 J±/2+l 



"fk = cos(fcao). 



The sums run over the wave vectors k from the lattice 
BriUouin zone —n/ao < k < n/aQ. 

Tio is the usual quadratic quasiparticle Hamiltonian of 
the linear spin-wave theory (LSWT) corrected by the di- 
agonal quadratic terms coming from a normal ordering 
of the quartic Dyson-Maleev bosonic interaction. It is 
important that these corrections (similar to Oguchi's cor- 
rections in antiferromagnetal3) renormalize the magnon 
excitation spectrum (and the ground-state energy) with- 
out changing its basic structure, i.e., the number of Gold- 
stone modes. In terms of quasiparticle operators the 
Hamiltonian Tin simply reads 
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FIG. 2. Ground-state energy of the (1, i) ladder as a func- 
tion of the interchain coupling J±. The solid line displays 
the first-order 0{1/S) spin- wave result for the ferrimagnetic 
state, Eq. (^. The dashed line indicates the ground-state 
energy of the disordered spin-singlet phase, as obtained from 
the LSWT. The Lanczos numerical results for a ladder with 
N = 10 rungs are denoted by open circles. At Jx = there is 
a well pronounced cusp in the data indicating a sharp transi- 
tion from the ferrimagnetic state into a disordered spin-singlet 
state. 



B. Perturbation series for A/a 



Cfe = ci(o- + l)(27fc + J±)/V^ -f- 2(2c27fc + C3 and 
D = 2ci(2 + Jx) + {a + l)(2c2 + JxCa)/^^. Here wf'''^ 
are the dressed quasiparticle dispersions. The functions 
i^^'^^ without 0(1/5) corrections will be referred to as 
bare dispersions. 

Finally, the quasiparticle interaction V includes two 
different terms: the quadratic bosonic interaction 



v. 



Y.^t^^if^i+^k~^^kPk, 

k 



(6) 



which is defined by the vertex functions 



rA±) Drik~Ck a-1 ( .-s 
T — —ci 1 7fe + — ) , (7) 



2£k 



and the quartic Dyson-Maleev interaction Vdm defined 
by the vertex functions V^^^ — ^12-34, ? = 1, . . . , 9.0 Here 
and in what follows we use the wave- vector abbreviations 
{ki,k2,k3,ki) EE (1,2,3,4). 



In spite of the fact that the expression for the ground- 
state energy (||) produces an excellent fit to the exact- 
diagonalization (ED) results in a large interval up to 
Jx = 10 (in Fig. we show a smaller interval), this 
fact by itself is not enough to make a statement about 
the accuracy of the SWT. It is more appropriate to use 
as a criterion the sublattice magnetization Ma, since 
the latter keeps information on the long-range spin cor- 
relations as well. The zeroth-order 0{V'^) result for 
niA = Ma/N simply reads = 5*1—01. We shall 
deal below only with the on-site magnetization mA, as 
the relation mA + ms = 81 — 82 {vfiB = Mb/N) is ful- 
filled up to an arbitrary order of the perturbation series. 
Note that in the antiferromagnetic case Si = ^2 the off- 
diagonal interaction V2 in Hb disappears and, as a result, 
there are no 0{V) corrections to rriA- In the ferrimag- 
netic case V2 7^ 0, so that the expansion for tua in powers 
of V has the form 



rriA = Si - ci - 



J_s;-VkVt^ + vt 

4SN ^ Ek f""" ■ 

k 



^k +^fc 



(8) 



In Fig. H we show results for ruA obtained from Eq. (||) 
by using subsequently the bare and dressed dispersions 
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OJ^^'^"^ from Eq. (|^). Comparing with the numerical re- 
sults, one indicates that the expansion in powers of 1/S 
(using bare dispersions) does not describe qualitatively 
the weak-coupling limit J± <^ 1: it predicts a small de- 
crease of TUA in the vicinity of J± — 0. On the other 
hand, we find that the expansion in powers of V (using 
dressed dispersions) gives a correct qualitative result in 
the weak-coupling limit J_l ^ 1 as well. The indicated 
problem of the standard 1/5 expansion might be a result 
of the enhanced fluctuations of the chain magnetizations 
about the common quantization axis, and/or of the usual 
zero-point fluctuations [as far as the extreme quantum 
system (1, ^) is concerned]. The above observation im- 
plies that at least in the extreme quantum case (1, ^) the 
expansion in powers of V is more reliable: such a view- 
point is strongly supported by the second-order 0{V^) 
result for jua (see the Appendix and Fig. ^. The above 
interpretation of the spin-wave series (as power series in 
V) is also adopted in the following analysis of the exci- 
tation spectrum. As a matter of fact, the difference be- 
tween the standard 1/S expansion and those in powers 
of V is mostly pronounced in the extreme quantum case 
(1, ^): for larger site spins both expansions practically 
coincide. 




FIG. 3. On-site magnetization (sublattice A) of the ferri- 
magnetic phase as a function of the interchain coupling J_l. 
The dashed and dashed-dotted lines display the series results 
up to first order in powers of 1/5 (bare dispersions) and V 
(dressed dispersions), respectively [Eq. (^]. The solid line 
shows the second-order 0{V^) series result for niA (see the 
Appendix). The Lanczos numerical results for the N = 12 
ladder are denoted by open circles. 

In the whole interval < < 3 shown in Fig. ^, the 
second-order result for tua differs by less than 0.3^ from 
the ED estimates based on the = 12 ladderE£l. For 
larger J± one clearly indicates a monotonic increase of 
the deviations from the numerical estimates, which can 



be ascribed to the enlarged reduction of the classical spin 
SSi — Si — mA'- in the limit Jj_ ^ oo the exact result is 



SSi 



It is shown in the next Section that the SWT 



description of the excitation spectrum exhibits similar 
precision but in a smaller parameter range {0 < J± < 2) 



III. EXCITATION SPECTRUM 



A. Free quasiparticles 



The quadratic Hamiltonian TCq defines two branches 
of spin-wave excitations {a and /3 magnons) described 
by the dispersion relations eI"''^^ = 2Slu'^'^\ Eq. (||), 
in the extended Brillouin zone — tt/oo < k < n/ao (see 
Fig. Ill) . The excited states composed of a and (3 magnons 
belong to subspaces with < {Si — S2)N — 1 and 
> {Si — S2)N + 1, respectively. The magnon ener- 
gies e\^'^^ are calculated up to the order 0{1/ S). In the 
long wavelength limit fcap 1, the ferromagnetic branch 
e\^^ has the quadratic Landau-Lifshits form 



_ Ps ,2 



0{k% 



(9) 



where Mq = {Si — S'2)/(2ao) and ps are, respectively, 
the magiie,tization density per chain and the spin-stiffness 
constantEj of the ferrimagnetic ladder. This form of the 
Goldstone modes is characteristic for Heisenberg ferro- 
magnets and reflects the fact that the order parameter 
(i.e., the feippmagnetic moment) is itself a constant of 
the motion.Ey An alternative approach, relying on the 
expected conformal invariance of the related XXZ ladder 
model, pify also be used to predict the quadratic form of 
Eq. (|).M 

The spin-stiffness constant ps and the magnetization 
density Mq play a basic role in the low-temperature ther- 
modynamics of ID models with a continuous symme- 
try and conserved order parameters.E3 Up to the order 
0(1/5"), Ps can be obtained from the Landau-Lifshitz re- 
lation d) and Eq. (|): 



ao5i52 



1 -ci- 



C2 



5(7 



s^' 



(10) 



At the zone boundary k — vr/ao the ferromagnetic 
branch exhibits an additional minimum so that in the 
vicinity of tt/oo the spectrum reads 



= At") + const (^^ - k 



(11) 



where Ai- '' is the excitation gap at the zone boundary. 
Such a branch folding of the excitation spectrum is typ- 
ical for ladder structures - itj-has been indicated in uni- 
form spin-i ladders as well.o The excitation mode at 
k ~ tt/oo reflects the underlying ladder structure and 
in the weak-coupling limit J± <ti 1 may be interpreted 
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as uniform rotations of the antiferromagnetic moment 



L = En=i(-1)"+' - S2(n)]. As ^ 0, the ex- 

citation gap A^"^ cx Jj^ goes to zero. For ferromagnetic 
couphngs J_i < 0, the k = tt/oo mode becomes unstable, 
thus producing a global instability of the ferrimagnetic 
state. The minimum at k = tt/uo persists up to the limit 

J± = oo: the related maximum of -E^,"'' changes its po- 
sition from k — 7r/(2ao) (for Jj^ = 0) to fc = vr/ao (for 
J± = oo). 

Similar folding appears in the optical antiferromag- 
netic branch Eif\ which can be characterized by the 



spectral gaps Ag''^ and A^-''^ at A; = and k = n, respec- 
tively. For instance, using Eq. (||), the spectral gap A^^^ 
up to 0(1/5") reads 



{2 + Ji_)iSi - S2) 1 



2c2 + c3Ji_ 
S^{2 + Ji 



(12) 



For the (1, i) ladder at Ji_ = 0, Eqs. (0) and (^ give 
the estimates Ps/{aQSiS2) — 0.761 and A, 



1.676. 



On the other hand, the LSWT Hamiltonian Hq pro- 
duces the parameters of the related classical system: 
Ps/ {aoSiS2) = 1 and Ag''^ = 1. The numerical estimate 
for the gap Ag''^ = 1.759 at Jj^ = 00 clearly demon- 
strates the importance of the 1/S corrections in Eq. (||). 




FIG. 4. Excitation spectrum of the (1, i) ladder for in- 
terchain couplings J± =0.1 and J± = 1. The dashed and 
sohd lines display, respectively, the series results up to the or- 
ders 0(y°) = 0(1/5) and C(y2). The 0(y°) results for the 
acoustic branch -E^"' (which are not displayed) closely follow 
the presented second-order curves. The symbols indicate our 
Lanczos numerical results. 



B. Quasiparticle interactions: second-order 



(«,/3) 



corrections for lj 



The quasiparticle interaction V produces corrections to 
the dispersions ^^1°'''^'^ only up from the order 0{V'^). The 
second-order corrections for the ferromagnetic branch 
tj^"' are described by the self-energy diagrams shown in 
Fig. 2 of Ref. |l3|. The respective analytic expression for 
the corrections reads 
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(13) 



Note that since the vertex functions Vj| , v/^-* 
and yili. 



(3) 

k ' ^ kp\pk ' kp\pk ' 

vanish at the zone center fc = 0, the 



gapless sliucture of 
diagramEj SW = 5{1 
function. 



(a) 



is preserved separately by each 
2 — 3 — 4) is the Kronecker 5 



The second-order corrections for 



come from simi- 



lar diagrams. The explicit expression reads 
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(14) 



The above expressions can be used to find the second- 
order corrections for the spin-stiffness constant ps and 
the excitation gaps aI"-*, Aq'^\ and Ai'^-'. The results for 
the ferrimagnetic state are summarized in Figs. and 
in Table | 
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FIG. 5. Gaps in the excitation spectrum of the (1, i) fer- 
rimagnetic ladder as functions of the interchain coupling J± . 
The solid lines represent perturbation series results up to sec- 
ond order in V. The open circles indicate Lanczos numerical 
results for a ladder with A'^ = 12 rungs. 



IV. ANALYSIS OF THE RESULTS AND 
CONCLUSIONS 

A comparison of the SWT and ED results (see Fig. |^) 
implies that already at the level of free quasiparticles the 
theoretical curves for the acoustic branch closely 
follow the numerical estimates almost in the whole Bril- 
louin zone apart from some vicinity of the wave vector 
k — 7r/(2ao). Note that the indicated discrepancy ex- 
ists only in the weak-coupling limit (the case J± — 0.1 
in Fig. Similar problem has been indicated for the 
(1, i) ferrimagnetic chainllj The point is that the second- 
order corrections for Eq. (|l3|), in the vicinity of 
k = 7r/(2ao) are very small (about 0.6% from the prin- 
cipal approximation for J± — 0.1). Therefore, in the 
weak-coupling limit Jj^ <C 1 the series does not describe 
in a proper way the acoustic branch close to the wave 
vector k = 7r/(2ao) (where the deviation is about 10% 
from the extrapolated ED result). Looking at the re- 
cently published Monte Carlo results for different {81,82) 
ferrimagnetic chainslla, it is seen that the discussed dis- 
crepancy in the region oi k — it / (2ao) is well pronounced 
only in the extreme quantum system (1, i). On the other 

hand, turning to the optical branch one finds an ex- 
cellent agreement between the second-order OiV^) per- 
turbation result and the numerical estimates (see Fig. ||). 
For instance, the theoretical result for the excitation gap 
l^f^ at = 0.1 differs by less than 0|^% from the ED 
estimate. Moreover, it can be expectedcj that thcr-third- 
order corrections further improve the above result .E3 The 



plots of the excitation gaps Ai-" , Aq , and Ajr (see 
Fig. H) can be used to find the range where the SWT 
still produces good quantitative results. The largest dis- 
crepancies are connected with the excitation g (at 
Jj_ = 1.5 the deviation from the ED result is 1.4%, but 
at = 3 it already exceeds 15%). The discrepancies 
grow with the interchain interaction, but up to Jx w 2 
the theory produces good quantitative results. 




0.0 0.4 0.8 1.2 1.6 2.0 
Jx 

FIG. 6. Spin-stiffness constants versus J± of different 
(Si, 5*2) ferrimagnetic ladders calculated up to second order 
of the perturbation series in V . 

The plots of ruA, and ps (see Figs. ^, ^, and 

also show that the short-range correlations, which are 
typical for the strong-coupling limit — > -l-oo, are qual- 
itatively established already at J_l =2. On the other 
hand, the same plots indicate a well pronounced crossover 
behavior for smaller interchain couplings (0 < J± < 1): 
the spin-stiffness constant ps and the on-site magneti- 
zation TUA reach their maximal values (respectively, at 
J± — 0.45 and 0.77), whereas the excitation gap Aq'^'' 
becomes minimal (at J± w 0.15). The discussed peculiar- 
ities appear as a result of the crossover between different 
types of short-range correlations without global change 
of the ground-state symmetry. 

Turning to the case of ferromagnetic interchain cou- 
plings {J± < 0), as mentioned above, one can expect that 
the model (|l|) is in a disordered spin-singlet phase. In the 
limit J± — > —00, the system is equivalent to the antifer- 
romagnetic spin-(5'i + 82) Heisenberg chain. Again the 
weak-coupling limit |Jj_| <C 1 is more interesting, espe- 
cially in the presence of exchange anisotropies. Here we 
shall restrict ourselves to short comments on the weak- 
coupling limit |Jj_| <C 1 of the isotropic (1,^) ladder. 
The ED results is Figs. ^ and || imply that already in the 
weak-coupling limit |J_l| <C 1 the short-range spin cor- 
relations qualitatively reproduce the corresponding spin 
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correlations in the spin-| antiferromagnetic Heisenberg 
chain. A quahtative LSWT analysis of the excitation 
spectrum predicts two degenerate branches of spin-wave 
excitations Ej^\ The lower branch E^^ ^ has a linear 
dispersion close to the wave vector k = 0, which is char- 
acterized by the spin- wave velocity 

^2)2+2(51-52)2 

On the other hand, -E^^"* describes gapped spin-wave ex- 
citations: the characteristic spectral gap at fc = is 

A = v^l J±P(5i + 52)2 + \J^\{Sf + 5|) + 4(5i - 52)2 . 

(16) 

It is seen that the relativistic modes remain stable (and 
the gap A is finite) up to the limit of noninteracting 
chains. The above LSWT result may be interpreted in 
the sense that the classical magnetic state in the region 



Jj^ < is swept out by the quantum fluctuations.EZI 
On the other hand, as is well known, a LSWT analy- 
sis cannot exclude the existence of gapped Haldane-type 
phases. To some extend, the discussed ferrimagnetic 
model {J± < 0) in the weak-coupling limit resembles the 
recently studied two-leg spin ladder constructed of ferro- 
magnetic spin--i cha ins w ith antiferromagnetic interchain 
couplings {J_i > 0).l3Ej It was found that for arbitrary 
J± > and isotropic exchange interactions the ground 
state is disordered and the system behaves like the spin-1 
Heisenberg antiferromagnetic chain. Although the above 
analysis cannot completely exclude the existence of a sim- 
ilar gapped phase for | Jl| ^ 1, we believe that the crit- 
ical phase remains stable up to J± = 0^. Such a conclu- 
sion is supported by the ED results which demonstrate 
well-established short-range correlations even in the limit 
I I <C 1 . In terms of the related nonlinear a model de- 
scription, on may expect that one and the same phase 
6 = IT (mod 27r) of the topological term will describe the 
whole parameter range J± < 0. 



TABLE L Results for the spin-stiffness constant ps, the on-site magnetization niA, and the excitation gaps (A 
and Ai''' ) of different J± — 1 ferrimagnetic ladders, as obtained from the second-order perturbation series in V . 

{Si,S2) Ps/{aoSiS2) rUA A^^^ A^' 

(1, 1) 0.8570 0.8484 1.4480 2.2303 3.6596 

(|,1) 0.8455 1.2517 2.9144 2.0687 4.9666 

(|,i) 0.9435 1.3992 1.3987 3.8191 5.2095 

(2,1) 0.9264 1.8285 2.8350 3.6952 6.5197 

(2, i) 0.9691 1.9236 1.3789 5.3623 6.7367 

(|,|) 0.9804 2.4382 1.3683 6.8885 8.2540 




FIG. 7. Short-range intrachain correlations in the 
spin-singlet phase of the (1, i) checkerboard ladder. The 
diamonds and circles represent, respectively, Lanczos nu- 
merical results {N — 10 ladder) for the spin-spin 
correlators ([Si(n) + S2(n)] ■ [Si(n + 2) + S2(n -I- 2)]> and 
([Si(n)-hS2(n)] ■ [Si(n + 4) + S2(ri + 4)]) The dashed lines 
indicate ED results (A'^ = 10) for the respective spin-spin cor- 
relators in the antiferromagnetic spin-| Heisenberg chain. 



In conclusion, the accomplished analysis clealy shows 
two basic implications of the underlying lattice structure 
for the ladder ferrimagnetic state: (i) the spin- wave exci- 
tations form folded acoustic and optical branches in the 
extended Brillouin zone and (ii) the ground state parame- 
ters (such as the on-site magnetizations niA and the spin- 
stiffness constant p^) show a crossover behavior in the 
weak-coupling region < J_i < 1. The presented results 
can be used in studies concerning the low-temperature 
thermodynamics of ferrimagnetic ladders. 
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FIG. 8. 



Short-range interchain correlations in the 



(1, 1) system, as obtained from the numerical diagonaliza- 
tion of the A'' = 10 ladder: (Si(n) • S2(n.))-filled circles, 
(Si (n) ■ 82(71 + 2))-open circles, (Si (n) ■S2{n + 4))-diamonds. 
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APPENDIX A: SECOND-ORDER 
CORRECTIONS FOR Ma 



The second-order correction for the on-site magnetiza- 
tion iriA can be expressed in the form 
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(Al) 

where (A)^^) denotes the second-order correction for 
the ground-state average of the operator A. The dia- 
grams connected to {alak}^^\ iPlPky^K and 
{akfik)'^'^^ are presented in Fig. ^. The related exphcit 
expressions read 
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FIG. 9. 



Second-order diagrams contributing to m^-*: 



(a), (b), and (c) diagrams give contributions, respectively, 
to (aIafe>^^^ (/3l/3fc and The diagrams for 

('^I/^I)^^^ can be obtained from the (c) diagrams by using 
the following vertex substitutions: c\: (7,2) (8,3), C2: 
(5,7) ^ (6,8), C3: (-,7) ^ (+,8), and C4: (+,4) ^ (-,4). 
The a (/3) magnon propagators are represented by solid 
(dashed) lines, whereas the vertex functions are denoted by 
their superscript indices. 
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